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NORMALIZING PARAMETERS FOR THE CRITICAL FLOW RATS 

OF SIMILE FLUIDS THROUGH NOZZLES 

by Robert C. Hendricks 

NASA Lewis Research Center 
Cleveland, Ohio 44135 U.S.A. 

SUMMARY 

Two-phase critical flow of simple fluids through a nozzle has been shown 
to nearly obey the principle of corresponding states. Quantum fluid depar- 
tures from the principle were resolved as a function of temperature for para- 
hydrogen and helium. 

The critical flow rates were normalized using 

G *- 

which nearly reduces the calculated equilibrium values of critical flow for 
classical fluids methane, oxygen, nitrogen and argon to a single isothermal 
curve. For the quantum fluids, p-hydrogen and helium, the normalizing param- 
eter becomes dependent on the isotherm and quantum fluid of interest. The 
normalizing parameter 

G [1 + ^q] 

where 4 'q = 0 for classical fluids, nearly reduces the critical flow rates 
for a! 1 simple fluids to a single isothermal curve. 

Experimental data for classical fluids nitrogen and methane for reduced 
pressures to 2.5 and reduced temperatures from 0.8 where taken and have been 
reported elsewhere. Experimental data for p-hydrogen, reduced pressures to 
4.5 and reduced temperatures from 0.87, were taken and analyzed herein. These 
data support the procedures for normalizing two-phase critical flow of simple 
fluids, including quantum fluids, through nozzles. 

INTRODUCTION 

The concept that there could be several thermodynamic processes which 
nearly obey the principle of corresponding states was advanced in Ref. 1 and 
extended to two phase choked flows through nozzles in Ref. 2. The critical 
flow rate for several simple fluids, c.f. methane, nitrogen, oxygen, and 
argon, were normalized using the parameter: 



’.-here the values of G* for several fluids are given in table I. The calcu- 
j-i'.-d 1 " normalized critical flow rate (G/^*) for the above four fluids were 
nearly collapsed to single isotherms. Experimental data for nitrogen and 


^Isentropic, equilibrium expansion of a fluid through a, nozzle was assumed; 
tnis does not advocate any particular theory but facilitates the establishment 
of normalization procedures. 
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methane verified the normalization procedure. 

The question of how to apply this procedure to normalize the critical flow 
rates of quantum fluids para -hydrogen and helium is addressed in this paper. 


SYMBOL LIST 




G* - 

P 
m 
T 
Z 
A 
e 


ft 


functions defined by 
eqs. 5 and 6 

critical mass flow 
rate 

normalizing param- 
eter for G 

pressure 

molecular mass 

temperature 

compressibility 

de Broglie wavelength 

intermolecular force 
constant 


density 


o 


empirical deviation function 
collision diameter 
mass flow rate 


Subscripts : 


c 

He 

o 

PHg 

P.G. 

R 

0,1 


ANALYSIS 


thermodynamic critical condition 
helium 

stagnation condition 
para-hydrogen 
perfect gas 
reduced 

reference conditions 


Gunn et.al. [3], found that quantum deviations in the reduced Joule- 
Thomson curve could be expressed in terms of the de Broglie wavelength 


A = 



( 2 ) 


The inversion curve was then expressed as the sum of the classical corre- 
sponding states part fQ (T R ) and a quantum deviation f^ (T R , A 2 ) 

Pr = f 0 ( t r) + f l < T R’ a2 > ( 3 ) 

In this report the physical thermodynamic critical constants, table I, 
will be used. This represents a departure from Ref. 3 where the constants 
were modified for quantum effects. 
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The isentropic equilibrium calculations for j>ar a -hydrogen and helium 
were performed, using the techniques of Ref. 4 and the thermophysical proper- 
ties procram GASP Ref. 5 for para-hydrogen and other fluids and HELP (inter- 
nal NASA-Lewis Research Center program) for helium. 


In Fig. 1, the reduced critical flow curve's for the classical fluids 
methane, nitrogen, oxygen and argon and the quantum fluids para-hydrogen and 
helium for Tr = 1.0 are given. Also presented are Tr = 0.8 and Tr = 1.2 
isotherms for nitrogen, para-hydrogen and helium. For the Tr = 1.0 isotherm, 
while the critical flow’ rate G varied nearly a factor of 10, the reduced 
curves for the above simple fluids deviate +_ 7 % at low pressure to l‘« at 
Pr = 4. This success implies that the normalization procedures of Ref. 2 can 
be extended to quantum fluids. Departures along other isotherms are pronounced 
and dealt with subsequently. 


Reduced critical flow rate plots for para-hydrogen and helium for selected 
isotherms were constructed, Figs.' 2 and 3, respectively. Fig. 4 represents 
the reduced critical flow rate plot for classical fluids (e. -. :.L,r •••.. an,, 
oxygen) . > 


Fig. 1 and overlays of Figs. 2 to 4 reveals considerable differences 
between para-hydrogen, helium, and the classical fluids except near Tr = 1. 
The problem now is how to resolve these differences. Assume that the reducing 
parameter for critical flow rate may be expressed as the sum of the classical 
reducing parameter and the quantum deviation: 


a * v 


( 4 ) 


It was found that while Eq. (4) is a function of both temperature and pressure, 
the pressure dependence appeared to be minimal see Fig. 1. Selecting a re- 
duced pressure of 3 (Pr = 3) as best representing the range of interest, one 
forms the following functions: 

(G/G*)pH 2 


f(I E ) = 

* - 1 

(G/G )classical 

(5) 

b(t e ) = 

( G/G*) He 

(6) 

1 

(G/G ) classical 


where classical refers to nitrogen, oxygen etc. Using Figs. 2 and 4, f(T R ) 
was found and plotted as Fig. 5. Similarly using Figs. 3 and 4, g(T R ) was 
found (see fig. 5). If the functions f(T R ) and g(TR) are valid, i.e. pres- 
sure independent, one should now be able to reduce the critical flow rate of 
say para-hydrogen and nitrogen to a set of single isotherms. To verify the 
I’uri'C L 'ii s f ( Tr ) arid r(Tp) i ho •Tit L*al fl vr f i para-Lyuv f . i, ana i, 
w<jre oal'-ulat.ed using 



reproducibility 


Of THE ORIGINAL P3B1 ,j| 
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Cp - 


u pH.. 

c 

U + f(Tp)]G* 


and Gp = 


J He 


[1 + g(Tp)lG* 


( 7 ) 


and the results are given as Fig. 6 along with -nitrogen. The a ;reement bctweer 
calculated isentropic equilibrium flow for classical fluids (dashed lines) 
and para -hydrogen ( -die elides) is quite good, with the exception of the region 
quite close to the saturation boundary for T_ < 1. Similar results were found 
for g(T R ), eq. (6) (crosses on Fig. 6). 

COMPARISON WITH DATA 


A limited number of tests were run with para-hydrogen in the blow down 
facility described in Ref. 4. The data and associated calculated parameters 
are given as table II. 

Comparing the normalized flow rate (G/G )pHg in the data section to that 
in the Calculated Parameters Section, reasonably good agreement can be found. 

A few selected data points from table II can be compared to the calculated 
isotherms of Fig. 2, see Fig. 7. As is the case with nitrogen and methane, 
the calculated isotherms lie above the data - up to 10# for Tp - 0.8, to 
being close at Tp » 1.; to 5-10# below for T R > 1 ; and good agreement pre- 
vails in the gaseous regime. 

In Ref. 2 critical flow rate data for nitrogen and methane are compared 
at selected isotherms. The equilibrium calculations serve as convenient 
reference lines and their use does not advocate any particular model; the 
comparison is between the data. Unfortunately the para-hydrogen data are not 
along equivalent isotherms, consequently no direct comparison with the data 
of Ref. 2 is made. The indirect comparison between calculations for classical 
fluids and experimental data, similar results for para-hydrogen and the agree- 
ment illustrated in Fig. 6, establish Eq. (7) as a technique to normalize the 
critical flow rate for simple fluids. 

CONCLUSIONS 

The normalization parameters for two phase flow of simple fluids through 
a nozzle has been investigated. Quantum fluid departures from the classical 
fluids were resolved as functions of temperature. 

The generalized normalization parameter for simple fluids is: 


where 



G 

[1 + i|tq]G* 



for classical fluids 
for para -hydrogen ' 
for helium 



'Mad* 





G 



P c , P c , Z c are the physical 
critical constants 


Along the isotherm = 1, 4^ is quite small and ~ G/G . While the 
critical flow rate G varied nearly a factor of 10, the reduced curves for 
simple fluids deviate + 7 % at low pressure to + 1$ at = 4. 

In general, the agreement between calculated characteristics and experi- 
mental characteristics for para-hydrogen is quite good. The prediction of the 
critical flow rate for nitrogen or oxygen and the prediction of critical flow 
rate for para-hydrogen also showed good agreement. 
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TABLE I. - CRITICAL CONSTANTS USED IN THE REDUCING PARAMETERS 


FLUID 

Pc 

Tc 

Pc 

„ _ Pc . 

G* 


MN/m 2 

K 

gm/cc 

“ c PcRT c 

z c 

gm/cm 2 sec 

Nitrogen 

3.417 

126.3 

.3105 

.2937 

6010.4 

Oxygen 

5.083 

154.78 

.4325 

.2922 

8673.9 

Methane 

4.627 

190.77 

.162 

.2889 

5093.7 

Argon 

4.865 

150.7 

.531 

.2921 

9404.2 

Para-hydrogen 

1.2925 

32.976 

.03143 

.3023 

1158 

Helium 

.22746 

5.2014 

.06964 

.3023 

724 


TA &X II . - HYDROGEN CHOKED FLOW DATA 


t c - 1.2926 m /»*, T . 31 » 976 Xf 
O c • .05143 pt/cc; /ePc 


CALCULATED PARAMETERS 


G 

g*/c*2»*c 


UfS 2*972 
1194 2.017 
1193 4*701 
U92 3.492 
USX 3.756 
1190 2*431 
1169 2*941 

1166 4,314 

1167 2*974 
1164 3*024 
1188 4.411 
1164 3,47 


.0294 
.0464 
.0337 
.0339 
*0907 ( 2.36 
.0361 1 1.62 
*0373 
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Figure 1. - Reduced critical »la* for several fluids along 
selected isotherms. 
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Figure 2. - Critical flow rate of para-hydrogen computed 
by Isentroplc equilibrium expansion. 
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Figure S. - Quantum deviation functions for critical flow 
of para-hydrogen and helium through a nozzle based on 
isentropic equilibrium expansion. 
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